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ABSTRACT 


Numerical calculations of the coherent and incoherent one-phonon scattering of slow neu- 
trons by a single KCl crystal are presented. The one-phonon scattering surfaces for a series of 
neutron energies and some different crystal settings have been determined on the basis of the 
eigen frequencies obtained by Tenerz when using a realistic model of KCl. Rough estimates of 
the seattering intensities have also been made. 


I. Introduction 


In 1938 and 1942 Pomeranchuk [1] and Seeger and Teller [2] discussed certain 
features of the inelastic scattering of thermal neutrons against crystals. Since then many 
authors have developed the theory in more detail. Especially two of them, Placzek 
and Van Hove [3], have emphasized that processes where only one vibration quantum 
participates would give valuable information on the crystal vibrations. In such cases 
measurements of the relation between the energies of the scattered neutrons and the 
scattering angle would give the frequencies for any single vibration wave vector, at 
least in principle. To make sure that multi-phonon processes do not disturb the meas- 
urements the energy of the incoming neutrons and the temperature must not be 
too high. In recent years some successful experiments have been performed to 
measure this energy exchange corresponding to one-phonon processes [4—10].1 

Squires [11] has made a calculation of the energy-angle relation for aluminium 
based on the Born-von Karman model for the crystal vibrations. The present paper 
contains a presentation of similar calculations for KCl at room temperature and a 
rough estimate of the intensity of the scattered neutrons is made. The eigen frequen- 
cies were determined by Tenerz [12] in connection with a study of the specific heat. 
The crystal was then considered a simple cubic one. In some respects these calcula- 
tions seem to be more reliable than a previous one by Iona, owing to the fact that 
the actual charge of the ions has been replaced by the effective charge, in agreement 
with experiments [13]. The short range forces were determined from the experimental 
data of the elastic constants and the residual ray frequency. Subsequently new cal- 
culations were made using a more detailed model and these results ought to be more 
trustworthy [14].? 

1 Dr. Larsson and collaborators at the Swedish Atomic Energy Company have also recently 
performed such experiments, with very good energy resolution. 

2 When the calculations presented in this paper were performed (1955) only the earlier results 
on the eigen frequencies were available. There are fairly large discrepancies between these and the 
later ones, not yet published, and therefore the numerical results in this paper must be taken only 
as a rather crude approximation for KCl. 
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II. The eigen vibrations of the crystal 


In the calculation of the frequencies mentioned above the K*- and Cl--ions were 
considered identical, apart from their charges, and a cubic reciprocal units cell with 
edge length 2/a was introduced, a being the distance between the nearest ions 
(a K+- and a Cl--ion). The frequencies and the polarization vectors of the vibrations 
were determined from a system of equations 


wo C,=> Diy C; (1) 
y 
where the dynamical tensor D,, depends on the wave vector q. For each q there are 
three different modes of vibration, characterized by 


w,(q) and C,(q), j=1, 2,3. 


Fig. 1. The direct unit cell; Fig. 2. The reciprocal unit cell. 
e K-atom, x Cl-atom. 


Two of these correspond to “‘transversal vibrations’? and one corresponds to a 
“longitudinal vibration’’.1 

The K and Cl nuclei scatter neutrons differently and therefore it would in our case 
be more convenient to consider the crystal as a face-centered cubic lattice (the unit 
cell is drawn in Fig. 1) occupied by a K+-ion at the lattice points and a Cl--ion in 
the middle of each unit cell. The corresponding reciprocal lattice is then a body- 
centered cubic one. Its unit cell, the first Brillouin zone, is drawn in Fig. 2, within 
the cubic reciprocal cell above.? 

The eigen vibrations are then determined from the following altered system of 
equations 


wuts. = > Dry Os kk 2D (2) 
k=1 refers to a K*-ion 
A Se ele let Bi Ds F553 


and normalized according to > |€*|?=1. For each q there are now six modes of 
k 


vibration, three acoustic ones and three optical ones, with their characteristic quant- 
ities 
O; (q), Cj (q), CF (q), = ib 2, eeey 6. 


L The polarization vectors are exactly parallel and perpendicular to the wave vector only for 
certain symmetry directions of q. 


* Note that a factor of 22 is included in our definition of the reciprocal lattice. 
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We note, however, that in this case the wave vectors fill a smaller volume than above 
(see Fig. 2). 

To use the numerical results by Tenerz we must identify the solutions from (1) and 
(2) respectively. Thus the three acoustic solutions of eqs. (2) correspond to the three 
solutions of eqs. (1) with the same wave vector and eigen frequencies. The relation 
between the polarization vectors is 


CF(q)=G(q)e"™*, j=1,2,3 (3’) 


where €,(q) refers to the representation in (1) and r, are the position vectors of the 
ions in the zero unit cell.t The optical branches in (2) on the other hand correspond 
to the solutions of (1) for a wave vector q’ outside the polyhedron. The frequencies 
and polarization vectors of the optical branches are determined from 


iq-r 
eum 


0; (q)=a,-(q'), CF (4) = ex Cy (q’) 3”) 
é=1, e=—1, 7=4,5,6, 7 =1,2,3 
where q' = q+ tT, t being a certain reciprocal lattice vector. 


III. The scattering surfaces 


The relation between the energy exchange and the scattering angle is determined 
from an energy-momentum law [1] (see Fig. 3). 


2m 
k? = ke ea; (9) 


k=k)+q-t 


Here k, = wave vector of the incident neutron. 
k =wave vector of the scattered neutron. 
<~ =reciprocal lattice vector, indicating Laue intensity maxima. 
€ =-+1, depending on whether neutrons lose or gain energy in the scattering. 
m =mass of the neutron. 


‘5 
Fig. 3. 
1 ry, = a(0,0,0), rp = a(1,1,1), where (1,0,0), (0,1,0) and (0,0,1) denote the unit vectors along the 
principal axis of the crystal. 
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Fig. 4a. Fig. 4b. 


Fig. 4c. Fig. 4d. 


Fig. 4a-f. Iso-frequency curves in the seattering plane. (The frequencies are arranged in such a 
way that the polarization vectors vary continuously.) 


a) Acoustic longitudinal branch. 
b) First acoustic transversal branch. 
c) Second acoustic transversal branch. 


d) Optical longitudinal branch. 
e) First optical transversal branch. 
f) Second optical transversal branch. 
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Via 


Fig. 4e. Fig. 4f. 


Our main task has been to solve (4) for a given Ky, thus determining k as a function 
of the scattering angle. The wave vectors of the scattered neutrons then define a 
certain scattering surface for each 7 and «, . 

In our calculations neutron energies corresponding to ky = 0, 1.8, 3.0 and 4.0- 108 
em-! (A=(22)/ky =00, 3.49, 2.09, 1.57 A) were considered and the scattering plane 
was supposed to be perpendicular to a principal axis of the crystal. Only those 
wave vectors q reaching points in the (xy0)- and (xy, +2/a)-planes of the reciprocal 
unit cell contribute to the scattering (see Figs. 2 and 3). 

The equations (4) were solved graphically in the following way.! On one sheet 
was drawn the contour of the reciprocal lattice in the scattering plane and into 
each cell were put series of iso-frequency curves, one series for each branch 7(@ = 
1, 2, ..., 41 1012 sec-!). Some of these iso-frequency curves are shown in Fig. 4. A 
blank cell indicates that the corresponding scattering has a vanishing intensity. On 
another sheet of transparent paper was drawn the vector ky, representing the wave 
vector of the incident neutron, and also a series of concentric circles with the radii 
k= Vie—e2 m/hw, where w has the same values as above. The latter sheet is super- 
imposed on the first one as shown in Fig. 4a, and the intersections between correspond- 
ing curves on the two sheets then determine the scattering surfaces in accordance 
with the energy-momentum law (Fig. 3). When a different direction of the incident 
neutron is required one need only rotate the k-sheet. 


IV. The scattering intensity 


The intensity of the scattered neutrons is obtained from the following expression 
for the differential scattering cross section (the cross section per unit solid angle 
and unit energy interval divided by # and per unit cell) [15]: 


1 The same method has also been suggested by A. Herpin, B. Jacrot, J. Phys. Radium, 16, 
35S (1955). 
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dot, Qa k x 
mh -— oly 
aidan ay Te, Dep & Ste? Sw? 
h (u- Cf )* (x- CF) coth (Bas) +e 
2NMao, 34 


- Whe- Wy pied : C= rp’) 


Op (x + q) 6 (@ — Wy + Es). (5) 


Here 2N =number of atoms in the crystal (V = number of unit cells). 
M =mass of the nuclei, taking the mean values for K and Cl. 


V =volume of a unit cell. 
hw, é 
p= akeT” where kz = Boltzmann’s constant. 
“Kp 


T' = temperature 


and s denotes the double index (q,7). 


he nk 
Se 2 Om yee 


x=k,—k 
<a,> are the coherent bound scattering lengths for K and Cl and their values are 
taken from Bacon [16]. 

<a> = 0.35: 10-12 cm, <aq)> = 0.99 - 10-12 cm. 


e “k are the Debye-Waller factors, which in the Debye approximation are determined 
from! 


@/T 
I? x” TK? x x 
2W,=2W,=2W= eco eae he 
Wa is aan Ue swice® (a) | (aa+5) (6) 
it 


and for © = 219°K, the Debye temperature of KCl, and 7’ = 293°K this simplifies to 
2W =2.45-10-18|k, —k|?. 
6(...) is an ordinary Dirac function and dp(x) = > (%—‘) means a periodic Dirac 


Tv 
function, singular at the reciprocal lattice points. 


If the polarization vectors are taken from (3’, 3’’) and normalized according to 
(2), the summation over k,k’ gives structure factors of the following type 


R= <a,> tite 


well known from the ordinary Bragg scattering. The value of this factor depends on 
whether the wave vector q is drawn from a corner lattice point or a centre lattice 
point of the body-centered reciprocal lattice and also on whether the eigen vibration 
considered is an acoustic or an optical one. Thus 


F=(axg>+aq>, |F|? =1.80- 10-24 cm? 
* See: R. W. James, The Optical Principles of the Diffraction of X-rays, p. 219, London 1948. 
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for an acoustic branch (7 = 1, 2, 3) and q belonging to a corner lattice point and also 
for an optical branch (j = 4, 5, 6) and q belonging to a centre lattice point. In all the 
other cases we have 


F=<ag)—<dq>, |}? = 0.41 -10-*4 em?. 
In (5) the summation over q is first transformed into an integration 


me 


Yay, 3 
q (22) 
(cell) 


. dq 


and this is immediately performed. The integration over w finally follows: 


= [o@ar- dm 


[ 0 @- e+ 60, (k-ky)) do = |-— 
dx 


where x = w — @) + €w,(K — Ko) 
da n-¢, 


= 1 
ai +é 


with n = unit vector along the scattering direction, 
v =velocity of the scattered neutron, 
e, = grad w,(q) = group velocity of the wave. 


A formula convenient for the calculations of the scattering intensity per unit solid 
angle and per cell is then obtained. 


si B aw) eho (ex C,)? coth (Ba) +e , (7) 
ky 2M on von, 2 


¢y is the unit vector along x and (q, x) are determined from eqs. (4). This formula 
was first derived by Waller and Fréman [17], and they emphasized the significance 
of the factor |1 + ¢n-e,/v|— for the scattering. 

In our calculations one of the acoustic and one of the optical branches have their 
polarization vectors perpendicular to the scattering plane, thus giving vanishing 
intensity (see Fig. 4). 

To obtain some numerical results on the intensity we must know the polarization 
directions and the group velocities. Unfortunately it was not possible to get more 
than a very rough estimate of the group velocities from the numerical material per- 
taining to the crystal frequencies. Nevertheless the calculated intensities, uncertain 
within perhaps a factor of two or more, give some valuable information on the 
scattering. 


V. The singular directions 


From (7) is seen that the intensity is infinite when the Laue scattering conditions 
are fulfilled (q =0, w, =0), as well as in the case of those scattering directions where 
(1+ ¢(n-e,)/v) vanishes. These singular directions have been studied earlier in more 


1 Note that according to the normalization of the polarization vectors cf we have | Cs |?=1. 
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detail under conditions where they occurred for small wave vectors (near Bragg 
scattering) [17, 18]. Here we will only show how to construct them when the scatter- 


ing surfaces are known. 
According to the energy-momentum law (4) we have 


dq=dk, dw=v-dk =—é grad w,(q):dq 
when k (and q) runs along a definite scattering surface. Thus 
(ec, +v):dk=0, ¢, =grad w,;(q) 


or in another form 
$¢;5b ViK=- vil; 


where n, is a unit vector normal to the scattering surface and y is a scalar factor. 
Multiplying both sides with 1/v n (n = the unit vector along the scattering direction) 
we get the factor appearing in (7) 


n:¢ 1 
1 = -n, 
+ £ P py (a N,) 


which vanishes only if k (k = kn) is a tangent of the scattering surface. The fact that 
y may be zero will not change this result. 


VI. The incoherent zero- and one-phonon scattering 


In addition to the coherent scattering there will be an incoherent one, due to spin 
dependence of the scattering lengths of the separate nuclei and to isotopic disorder 
in the crystal. As for the coherent part the incoherent scattering may be divided into 
zero-phonon, one-phonon etc. scattering, according to the number of vibrational 
quanta involved. Only the zero-phonon part gives discrete energies of the scattered 
neutrons (elastic scattering), whereas in the one-phonon processes the scattered 
neutrons have for all angles a continuous energy spectrum extending from h(a, — wn) 
to h(w» + w,), and analogous for the higher phonon processes. w,, means the maximum 
vibrational frequency. In principle it would be possible to distinguish experimentally 
neutrons scattered in an incoherent zero-phonon and a coherent one-phonon process 
from those scattered in other processes, using the fact that in the first case the 
neutrons have discrete energies. However, the detectors have a finite resolving power 
and this limits the practical use of this distinction. It would therefore in some cases 
be necessary to correct for these other processes. 

The differential scattering cross sections per unit cell are written below for the 
incoherent zero- and one-phonon processes [19].1 The two-phonon scattering, espec- 
ially the coherent one, is discussed in a separate paper [20]. 


(8) 


1 See also Placzek and Van Hove, loc. cit. 


206 


ARKIV FOR FysIkK. Bd 13 nr 17 


d E 
0.5 

vo a a . 
mie cc Piggleaia 
Ser antansfantfanti 

jee i 
-10 —0.5 0 0.5 10 
g (®) 


Fig. 5. To get the normalization according to formula (13) the values obtained from the figure 
must be multiplied by the factor 0.976. 


e *“ is, as above, the Debye-Waller factor and A means the sum of the incoherent 
differential scattering cross section of K and Cl taken for bound nuclei [16]. 


A = 0.26- 10-74 cm? 
g(@) is defined by 


g(w)=5 at [coth (Bw)—1], |w|<1 


=0, |o|>1 (10) 
where /(w) denotes the vibrational frequency function with the neutron “‘frequency”’ 


inserted instead of the vibrational one. For negative frequencies f(w) is defined as 
follows 


f(—@) = f(a) . (11) 


and normalized to unity! 
1 
[f()do=1 (11) 
y(0) is related to the Debye-Waller function by 
he 
2W =o?) (12) 


1 ~ is here measured in units of the maximum vibrational frequency w,, = 4.1 + 10'* secq?. 


207 


A. SJOLANDER, On the scattering of slow neutrons by KCl 


and is introduced to normalize g(w) as follows 


1 


[g(o)dw=1. (13) 


-1 


On very general grounds we may conclude that f(w) is everywhere continuous 
whereas df/d@ has a certain number of discontinuities [21, 22, 23], the strongest 
ones being of the following type: 


f(w)=|m@—a,|?, w>a, 
=0, W< We (14) 


g(@) is shown in Fig. 5. 


VII. Results on the scattering surfaces 


The numerical results on the one-phonon scattering surfaces are presented in a 
series of figures. The direction of the incident neutrons is first along a principal 
axis of the crystal, and in Fig. 10 some results for different incident directions are 
given. By using full and dotted lines it has been possible to give an indication of the 
variation of the intensities in the same figures, dividing the strength into three 
classes: 

0 — 3 10-* barns/steradian. cell: weak intensity 
3 — 15 10-? barns/steradian. cell: medium intensity 


> 15 10-? barns/steradian. cell: strong intensity. 


The concentric circles, indicating the lengths of the neutron wave vectors, are 
drawn to make a reading off easier and they are labelled both by the corresponding 
vibrational frequency w, measured in units of 102 sec-1, and by the corresponding 
neutron energy h*k*/2m. It must be noted that the intensities are very approximate. 
However, they give an idea of the variation with the scattering angle. 

Fig. 6, kg=4.0-108 em™’ (4=1.57 A): Fig. (6a) shows the scattering lines and 
Fig. (6b) shows the intensities in more detail for the most intense lines. To indicate 
which curve in (b) belongs to a certain line in (a), the corresponding curves in the 
two figures are labelled by the same letter. The strong elastic peak corresponds to 
the incoherent zero-phonon scattering. The Laue scattering conditions are nearly 
fulfilled in certain directions and in (b) we see how the intensity increases when 
approaching a Bragg scattering angle. 

Fig. 7, ky=3.0-108 cm~* (A=2.09 A): Especially the scattering lines representing 
phonon emission have been changed considerably compared with those in Fig. 6. 
Closed curves have been developed and by drawing the tangents of these from the 
footpoint of the incident wave vector we can construct the singular points defined 
by the vanishing of |1+¢(n-e,)/v|. 

Fig. 8, kj =1.8-108 em=! (A=3.49 A): No indication of the intensity is given. Note 
how the scattering surfaces for phonon emission disappear when the neutron energy 
decreases. At the same time the remaining surfaces are more separated from each 


Hon and the variation of the neutron energy with the scattering angle is more 
obvious. 


208 


ARKIV FOR FYSIK. Bd 13 nr 17 


6 — strong int! 
ee 14 ---medium oA 

22 eseree weak Jur 
0.053 30f 

38 


Fig. 6. ky = 4.0- 108 em™ (A = 1.57 A). 
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Fig. 7. ky = 3.0- 108 em- (A = 2.09 A). 


Fig. 9, kj =0 (A= co): For cold neutrons the wave length is large and then the 
scattering surfaces are approximately the same as those for infinite wave length. 
The intensity tends to infinity as 1/kp. 

Fig. 10, ky=4.0-108 em™* (A= 1.57 A): The direction of ky is varied in a principal 
plane, and it is shown how the scattering surfaces are altered within a certain 
reciprocal cell when passing a lattice point. To make Fig. (a) clearer the lines 
corresponding to w, and w, are drawn in separate cells. The angle between a principal 
axis and ky is indicated and thus the lines for 0° correspond to those in Fig. 6a. In 
(b) more detailed results are given for the w,-lines, including an indication of the in- 
tensity. The factor |1+¢ (n-e,)/v| has then been excluded and therefore the real 
intensity is appreciably increased when approaching a singular point. 
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